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SENSITIVITY OF MAGNETOSPHERIC MULTI-SCALE (MMS) 
MISSION NAVIGATION ACCURACY TO MAJOR ERROR 
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Corwin Olson,* Anne Long,* and J. Russell Carpenter* 


The Magnetospheric Multiscale (MMS) mission consists of four satellites flying 
in formation in highly elliptical orbits about the Earth, with a primary objective 
of studying magnetic reconnection. The baseline navigation concept is 
independent estimation of each spacecraft state using GPS pseudorange 
measurements referenced to an Ultra Stable Oscillator (USO) with 
accelerometer measurements included during maneuvers. MMS state estimation 
is performed onboard each spacecraft using the Goddard Enhanced Onboard 
Navigation System (GEONS), which is embedded in the Navigator GPS 
receiver. This paper describes the sensitivity of MMS navigation performance to 
two major error sources: USO clock errors and thrust acceleration knowledge 
errors. 

INTRODUCTION 

The Magnetospheric Multiscale (MMS) mission consists of four satellites with identical 
instrument suites flying in formation in highly elliptical orbits about the Earth. The primary 
objective of the mission is to study the phenomenon known as magnetic reconnection, which is a 
process that converts magnetic energy into heat and kinetic energy of charged particles. After 
launching into highly elliptical orbits and perigee raising maneuvers to place the spacecraft in 
orbits with a 1.2 Earth Radii (Re) perigee and 12 Re apogee, the mission begins with the 
commissioning phase in which initial testing of all spacecraft systems will occur. Phase la 
begins when the orbit has evolved such that the apogee is on the day side of the Earth. In Phase 1 
the formation separations will be adjusted from 160km to 10km in an effort to study reconnection 
on the dayside boundary between the solar wind and the Earth’s magnetosphere known as the 
“bow shock” region. The following period when the apogee is once again on the night side of the 
Earth is Phase lx, after which Phase lb starts when the apogee has returned to the dayside. 
During Phase 2a, a series of apogee raising maneuvers are executed to incrementally increase the 
orbital apogee of each spacecraft from 12 Re to 25 Re. Phase 2b will begin when the apogee 
raising for all four spacecraft is complete. During Phase 2b the formation separations will vary 
from 400km to 30km, with the optimum separation chosen by mission scientists after the full 
range of different formation sizes has been executed. During Phase 2b the mission will study the 
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night-side reconnection events that occur in the Earth’s magneto-tail. After Phase 2b, the mission 
will end and the spacecraft will be de-commissioned. Reference 1 contains a more detailed 
description of the mission concept. 1 

Figure 1 shows the high-level navigation concept for MMS along with the configuration of 
each spacecraft’s navigation and attitude control elements. Each spacecraft’s state is 
independently estimated using Global Positioning System (GPS) pseudorange measurements 
referenced to an Ultra Stable Oscillator (USO). MMS state estimation is performed onboard each 
spacecraft using the Goddard Enhanced Onboard Navigation System (GEONS), which is 
embedded in the Navigator GPS receiver. The Navigator GPS receiver employs a new weak 
signal tracking technology that significantly improves reception of GPS signals when the 
spacecraft is above 10 Re. 2 Thrust acceleration measurements from an accelerometer within the 
Attitude Control System are included during maneuvers. The estimated state consists of the 
position, velocity, clock bias, clock bias rate, and clock bias acceleration. The estimated states 
are periodically downlinked and used in the MMS Flight Dynamics Operations Area (FDOA) to 
generate definitive and predictive products to support mission and science operations. 
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Figure 1 Navigation Concept and Onboard Configuration. 

Navigation knowledge error, or the difference between the estimated state and the truth state, 
can be attributed to a number of different sources. The primary contributors include the 
truncation of the onboard gravity model to 8x8 due to onboard processing limitations, solar 
radiation pressure (SRP) modeling errors, GPS pseudorange measurement noise, GPS ephemeris 
errors, receiver clock errors, and thrust acceleration knowledge errors. This paper examines the 
sensitivity of MMS navigation performance to USO clock errors and thrust acceleration 
knowledge errors. An understanding of this sensitivity is needed to ensure that all mission 
requirements are met, especially if the USO clock or accelerometers do not perform as expected 
during the mission. 
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The MMS mission imposes several challenging navigation requirements due to the fact that 
the majority of the orbital period occurs above the GPS constellation where GPS signal 
acquisition is sparse. Critical requirements include the following 3 : 

• A maximum USO clock bias error of 325 microseconds is needed to meet the science 
requirement for a maximum relative clock bias error of 1 millisecond between 
spacecraft. 

• A maximum definitive Brouwer-Lyddane mean semi-major axis (SMA) error of 100m 
above 3 Re is needed to meet formation maintenance maneuver planning 
requirements. 

• A maximum relative definitive position error of 1% of the separation or 100m, 
whichever is larger, is needed to meet science requirements. 

• A maximum relative predicted position error growth of 200 m/day is needed to meet 
collision avoidance requirements. 

The USO clock error model parameters that were varied in this analysis include the clock 
process noise parameters and the amplitudes of the aging affect and umbra-induced temperature 
variations. The thrust acceleration knowledge error model parameters that were varied in this 
analysis include the direction bias, magnitude bias, magnitude drift due to random walk, 
magnitude scale factor, and temperature variation. These values are varied in an effort to 
determine the sensitivity of absolute and relative position and velocity navigation error and clock 
bias errors to these error sources. 

This analysis uses as a baseline a segment of the mission Phase 2b. Phase 2b is the most 
challenging phase of the mission for navigation because there are significantly fewer GPS 
measurements acquired for a much greater percentage of the orbit, and the receiver must 
propagate through over 2.1 days of the 2.8 day orbit without GPS measurements during apogee 
due to onboard power constraints. During the simulation time span, the nominal formation scale 
size (approximate distance between spacecraft during the science region of interest) is 30 
kilometers. The simulation includes a pair of formation maintenance maneuvers, which are 
designed to maintain the quality of the spacecraft formation. The first of these two maneuvers 
occurs immediately after exiting the science Region-of-Interest (ROI), which is defined as the 
region above 1 5 Re for Phase 2b. The second maneuver is immediately before re-entering the 
science ROI, after traveling through perigee. A pair of Phase 2b formation maintenance 
maneuvers for a 400km formation size scenario along with the final two maneuvers of the apogee 
raising sequence in Phase 2a are also examined to determine the sensitivity of navigation error to 
acceleration errors during a longer maneuver. Due to page length limitations, the Phase 2b 
400km separation and Phase 2a results are only briefly summarized below. 

SIMULATION METHODOLOGY 

Figure 2 outlines the simulation methodology that was employed for this analysis. The 
Mission Design (MD) code is used to generate the truth trajectory and the truth finite maneuvers, 
and the truth trajectory is fed into the Measurement Data Simualtion (DatSim) program to 
generate GPS pseudorange (PR) measurements. The accelerations from the truth finite 
maneuvers are modified to include knowledge error. The measurements and accelerations with 
knowledge error are ingested into GEONS, which generates estimates of the state and covariance 
values for the state. These states and covariance values are used along with the truth ephemeris to 
compute errors and verify requirements. 
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Figure 2. Simulation Methodology. 


TRUTH TRAJECTORY GENERATION PROCEDURE 

Figure 3 illustrates the procedure used to generate the truth trajectory and maneuver accelerations. 
The End-to-End (ETE) code referenced in the procedure is a code set designed to simulate all 
phases of the MMS mission from commissioning to the end of Phase 2b as one continuous 
simulation in order to verify all mission design and Conjunction Assessment (CA) requirements. 



Figure 3. Truth Trajectory Generation Procedure. 


TRAJECTORY MODELS, MEASUREMENT MODEL, FILTER SETTINGS 

Table 1 lists the primary models and associated values used to propagate the truth and filtered 
trajectory of each MMS spacecraft. The truth trajectory was generated using the FreeFlyer 
mission analysis software application. The GEONS filter models are consistent with the models 
being used to support the MMS mission, which are tailored to execute within the limited 
resources of the flight processor within the Navigator receiver. 
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Table 1. Trajectory Propagation Models. 


Simulation 

Parameter 

FreeFlyer Truth 

GEONS Filter 

Nonspherical Earth 
Gravity Model 

21x21 Joint Gravity Model-2 

8x8 Joint Gravity Model-2 

Point Mass Gravity 

Sun, Moon using DE 405 ephemeris 

Sun, Moon using analytical fit to 
DE 406 ephemeris 

Atmospheric Drag 

Analytical constant temperature 
hydrostatic model, C D of 2.2, Drag area 
of 7.1 m 2 

Analytical fit to Harris Priester model, 
C D of 2.2, Drag area of 7. 1 m 2 

Solar Radiation 
Pressure 

Spherical model, C R of 1.8, SRP area of 
2.026712 m 2 

Spherical model, C R of 1.8, SRP area of 
2.02 m 2 

Integrator 

8(9) Variable Step Runge-Kutta 

4 th Order Fixed Step Runge-Kutta 

Integration Stepsize 

10 seconds 

30 seconds 

Maneuver Model 

Finite burns 

Accelerometer measurements averaged 
over 10 seconds, including acceleration 
knowledge errors 


Table 2 lists the simulation parameter values used for measurement simulation. These values 
were selected to match the measured performance of the Navigator GPS receiver for specific 
ranges of the signal-to-noise-ratio (C/NO) value. 


Table 2. Measurement Acquisition and Error Models. 


Simulation Parameter 

Nominal Values 

Measurement Type 

GPS Pseudoranges (PR) 

Navigator Receiving Antenna 

Composite torroidal model for spinning spacecraft with peak gain 
normal to the spin axis 

GPS PR Measurement Rate 

1 Measurement set every 30 seconds for each formation member, 
with measurements from a maximum of 12 visible GPS SVs when 
below 15Re, no measurements above 15Re in Phase 2b 

GPS Acquisition Threshold 
(Based on Navigator 
Performance) 

C/N0 > 40 dB-Hz: 95% probability of acquisition, minimum 
acquisition delay of 600 sec 

C/N0 = 28 to 40 dB-Hz: 75% probability of acquisition, minimum 

acquisition delay of 600 sec 

C/N0 < 28 dB-Hz: 0% probability of acquisition 

GPS Antenna Model 

GPS IIA with Transmitting Antenna Effective Isotropic Radiated 
Power (SVEIRP) =13.4 dB-Hz 

GPS Constellation 

21+3 GPS Space Vehicles 

GPS Space Vehicle (SV) 
Selection 

Min to max Transmitter ID values 

la GPS PR Noise (Based on 
Navigator Specifications) 

4 meters for C/N0 > 40 dB-Hz (strong signal) 
10 meters for C/N0 < 40 dB-Hz (weak signal) 

GPS Ephemeris and Clock 
Errors 

Sinusoidal ephemeris model with 2 meter amplitudes for each 
component 

Ionospheric Delay Model 

GPS Ionospheric Model based on GPS broadcast coefficients 


Table 3 lists the parameters incorporated in the GEONS filter, which is a UD factorized 
extended Kalman filter. 4 
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Table 3. GEONS Filter Settings. 


Parameter 

Nominal Values 

Estimation State 

Position, velocity, clock bias, clock bias rate, 
clock bias acceleration 

Initial position and velocity state errors 

300 m and 15 mm/sec per axis (la) 

Initial Clock bias, drift, and acceleration errors 

150 km (500 jus) (la), 46 mm/s, 2.3 pm/s 2 

GPS PR standard deviation 

40 m 

Minimum Height of Ray Path Altitude 

1000 km (eliminates measurements with 
largest ionospheric delays) 


SENSITIVITY TO USO CLOCK ERRORS 

The following analysis was conducted to determine the sensitivity of onboard clock bias, drift, 
and acceleration estimation for the MMS mission to various clock stability parameters. These 
parameters include clock process noise parameters and the amplitudes of the aging effect and 
Umbra Temperature Events. The effect of aging refers to the change in the clock bias 
acceleration as time passes after the USO is turned on. Umbra Temperature Events are sinusoidal 
variations in the clock bias acceleration that have a period equal to the orbital period and are 
applied for every orbit to represent the worst case where for every orbit the spacecraft enters the 
umbra of the Earth. An understanding of this sensitivity is important due to the 325 microsecond 
requirement on clock bias knowledge. This requirement is most challenging during periods 
without GPS measurements when the clock bias must be propagated. 

These parameters are set by the user in the DatSim program, which generates measurements 
for the simulation. 5 Process Noise Variance Rates (PNVRs) for the clock parameters were 
modified as necessary in GEONS to ensure that the filter used to estimate the state was properly 
tuned. The USO clock error model is described in detail below. 


USO Truth Clock Model 


The truth clock model incorporates the following parameters: initial clock bias, initial clock 
drift, initial clock acceleration, frequency white noise, frequency random walk, frequency random 
run, aging model amplitude and time constant, and the amplitude, start, and stop times for 
periodic variations due to umbra-related temperature changes. 

The clock bias acceleration due to the aging effect is described in the equation 


U aging 


(',) = 




( 1 ) 


where the coefficients b x and b 2 are determined by fitting to empirical data to characterize the 
magnitude and effective rate of decay, t { is the current time and t 0 is the initial time. 

Solar and lunar umbras produce temperature variations during the orbital period in which the 
umbra event occurs. The effect of frequency variations that occur due to temperature variations 
before, during, and after umbras are modeled by including a clock bias acceleration term for each 
umbra interval. This term is modeled as a sinusoidal variation, with amplitude equal to b 3 and 
period equal to the duration of the event: 
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where f x is the change in the fractional frequency per 1°C change in temperature, and P is the 
period that is equal to the duration of the event, modeled as one orbital period. 

The full truth clock propagation is performed sequentially as follows: 
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where AT is the clock prediction step size in seconds and Q is the clock noise matrix, defined 
below. The values s 1 , s 2t and s 3 are normally distributed deviates with zero mean and unit 
variance, which can be initialized with the same or different random number seeds. 

For the correlated clock error model used in this analysis, the components of the clock noise 
matrix Q is defined such that QQ T is equal to the full clock covariance matrix: 
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and where q lt q 2r and q 3 are process noise variance rates associated with frequency white noise, 
frequency random walk, and frequency random run, respectively. The input simulation 
parameters h 0 and h_ 2 are white noise and random walk Allan variances that are related to the 
process noise variance rates q 2 and q 2 as shown in the equations 



q 2 = 2n 2 h 


The baseline scenario for this analysis was the MMS Phase 2b, 30km separation case. Table 4 
lists the baseline clock error parameters used in the measurement simulation software for each 
satellite for Phase 2b. These values are consistent with the expected performance of the USOs 
being procured for MMS. Table 4 also shows the clock model values used in the GEONS filter 
for this study. 

Table 4. Baseline Clock Error Model Parameters. 


Simulation Parameter 

Truth Model Value 

GEONS Filter Value 

Initial Bias 

3.6343x10° s 

3.0655x10° s 

Initial Rate 

1.54xl0' lu s/s 

0.0 s/s 

Initial Acceleration 

2.0x1 O’ 15 s/s 2 

0.0 s/s 2 | 

Frequency White Noise qi 

1.2E10' 22 s'/s (ho=2.4xlO' zz ) 

Not Modeled 

Frequency Random Walk q 2 

1.58x10-“ s'/s 8 (h. 2 =8xl0' is ) 

Frequency Random Run q 3 

10° 8 s'/s 5 

Aging Amplitude b ls 
Time Constant b 2 

5x1 O' 9 , 10 days (from turn-on) 

Periodic Variation Amplitude 
b 3 = (27i/P)10' n (AC) (s/s 2 ) 

3.7 h umbra in Phase 2b: 2.0 xlO' 15 


Variations Performed 

The clock parameters were varied in a systematic method to examine the sensitivity of 
onboard position, velocity, clock bias, clock drift, and clock acceleration estimation performance 
to the various clock stability parameters. The parameters varied include the noise parameters h 0 , 
h_ 2 , and q 3 , as well as the amplitude of the aging affect b l and the Umbra Temperature Events b 3 . 
Process Noise Variance Rates (PNVRs) for the clock parameters were modified as necessary in 
GEONS to ensure proper tuning. 

Table 5 below contains the variations of the clock model that were conducted for this analysis. 
The first run was with the baseline values that were consistent with the expected performance of 
the USOs being procured for the MMS Navigator receivers. Next each term was increased in 
magnitude to assess the sensitivity to degraded performance. 

Table 5: Baseline and Successive Clock Model Variations 


Case 

ho (s'/s) 

h .2 (s z /s J ) 

q 3 (sVs>) 

b, (s/s) 

b 3 (s/s') 

Baseline MMS clock 

2.4E-22 

8E-28 

IE-38 

5E-09 

2E-15 

o 

o 

* 

o 

43 

2.4E-20 

8E-28 

IE-38 

5E-09 

2E-15 

h. 2 * 100 

2.4E-22 

8E-26 

IE-38 

5E-09 

2E-15 

q 3 * 100 

2.4E-22 

8E-28 

IE-36 

5E-09 

2E-15 

b, *2 

2.4E-22 

8E-28 

IE-38 

10E-09 

2E-15 

b 3 * 2 

2.4E-22 

8E-28 

IE-38 

5E-09 

4E-15 

All values increased 

2.4E-20 

8E-26 

IE-36 

5E-09 

4E-15 
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Phase 2b Orbit Radius and GPS Visibility 

During Phase 2b, the MMS spacecraft receives far fewer GPS measurements when it is above 
the GPS constellation than when it is below the GPS constellation and no measurements at all 
when it is above 15 Re due to onboard power constraints. This lack of measurements results in 
propagation over long periods of time which can lead to large state errors, especially if maneuvers 
have occurred. As a result, Phase 2b is the most challenging phase of the mission for navigation 
of the spacecraft. 

Figures 4 and 5 display the MMS satellite radial distance from the Earth and the number of 
GPS satellites (SVs) visible during the simulated segment of Phase 2b. These figures show the 
large variation in the number of visible GPS satellites as a function of satellite radial distance. 
When compared with plots of estimation error results shown in the next section, the relationship 
between the number of GPS satellites tracked by the receiver and the estimation state error is 
apparent. 
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Figure 5. Number of GPS Satellites Visible for Phase 2b. 


Clock Model Sensitivity Results 

All results below are for spacecraft MMS 1 , which is representative of the sensitivity trends 
seen for the other spacecraft. The first case was run using the baseline clock model as described 
in Table 5, the results of which were compared to each variation described in Table 5. Figure 6 
shows the clock bias, drift, and acceleration error and formal standard deviation estimated by the 
filter with baseline clock settings as defined in Table 5. Comparing Figure 6 to the GPS visibility 
plot in Figure 5, the clock state error starts to grow when GPS measurements are no longer 
available and continues to grow until measurements are once again received, at which point the 
error drops dramatically to much smaller levels. The periodic variations in the clock acceleration 
error are due to the presence of umbra variations in the clock truth, but not in the clock estimation 
model. 
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Figure 6. Clock Bias, Drift, Acceleration Error and Covariance with Baseline Clock Settings. 

Figures 7 and 8 show the RSS position and velocity errors along with the filter’s formal 
standard deviation for the baseline case defined in Table 5. The points labeled Ml and M2 are 
the times of the first and second formation maintenance maneuvers. Notice that the error 
increases significantly immediately after M2 because after knowledge errors are induced by M2 
the estimated state must be propagated for about 50 hours through apogee until GPS 
measurements are once again acquired by the receiver. 
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Figure 8. Velocity Error and Covariance with Baseline Clock Settings. 

It was determined that, of the initial variations, only the cases that increased h_ 2 and b 3 
produced significant changes in the clock bias, drift, and acceleration error. All cases met the 325 
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microsecond clock prediction accuracy requirement, but if the h. 2 parameter is increased by 
another order of magnitude to 8xl0' 25 , the requirement is violated. This indicates that the h_ 2 
parameter has the smallest allowable variation in order to meet the requirement. 


Figure 9 shows the clock bias, drift, and acceleration error and covariance with the Random 
Walk Allan Variance h_ 2 clock parameter increase as defined as h_? * 100 in Table 5. It is clear 
that the increase of the h_ 2 parameter introduces significantly more bias, drift, and acceleration. 





Figure 9. Clock Bias, Drift, Acceleration Error and Covariance with h. 2 Increased to 8x10 26 . 

Figures 10 and 1 1 show the RSS position and velocity errors along with covariance for the h_ 2 
clock parameter increase as defined in Table 5. When the filter is appropriately tuned, this case 
did not have a significant impact on the position and velocity errors over the propagation period 
immediately following the second maneuver, as the additional error introduced was absorbed by 
the clock states. 
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Figures 12, 13 and 14 below illustrate the maximum clock bias, RSS position, and RSS 
velocity errors that result from the nominal baseline case and each variation. For all cases, the 
maximum position and velocity errors occur after the second maneuver. It is clear that the h. 2 
increase case, as well as the b 3 increase case, produce significantly larger levels of clock bias 

error as compared to the baseline. Position and velocity errors are not sensitive to these changes. 
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Figure 14. Maximum RSS Velocity Error for Each Primary Clock Error Parameter Variation. 


SENSITIVITY TO THRUST ACCELERATION KNOWLEDGE ERRORS 

The MMS onboard accelerometers will generate 10-second averaged acceleration 
measurements that GEONS will process during each maneuver, hi order to simulate these 
measurements, an acceleration knowledge error model was developed for each formation 
maintenance and apogee raising maneuver that includes errors that are consistent with the 
expected “worst case” acceleration knowledge errors. 

Thrust Acceleration Knowledge Error Model 

The Guidance, Navigation and Control (GN&C) subsystem will convert the raw acceleration 
measurements in the body frame to the mean of J2000 inertial (MJ2000.0) frame using the 
spacecraft attitude, and then average these values over 10 seconds to the nearest whole 10 second 
mark. Table 6 below describes the thrust acceleration knowledge error model parameters, which 
were derived from analysis conducted by the MMS GN&C subsystem team, and which include a 
50% modeling uncertainty factor. 
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Table 6. Thrust Acceleration Knowledge Error Model Parameters. 


Parameter 

1<t Value 

Direction Bias [dir mag ) - 
Formation Maintenance Maneuvers 

0.681° 

Direction Bias [dir mag ) - Apogee 
Raising Maneuvers 

0.447° 

Magnitude Bias ( bias mag ) 

1-5 Hg 

Magnitude Drift [drift mag ) 

22.5 pm/s/s 1/2 

Magnitude Scale Factor 

[ScaleFactor) 

0.01625% 


The thrust acceleration knowledge error model provides the perturbed acceleration vector 
expressed in the MJ2000.0 frame averaged over a 10 second interval, a mth eiror , given by 

a Klt h_e,ro, = Upertmertml (l^ ||0 ^ SeC ) + A Ve f T mag )/(10sec) (9) 

where u pert inertia] is the direction of the perturbed acceleration vector in the MJ2000.0 frame, ||5|| 
is the magnitude of the true acceleration, and A Verr mag is the total delta-v magnitude error, as 
detailed below. 

The total delta-v magnitude error is computed as the sum of the following contributions: 

A Verr mag (t) = A Verr magbias + AVerr magScaleFactor {t) + AVerr mag TempVar (t) + AVerr drift (t) (10) 

The delta-v magnitude bias (a Verr mag bias ) which will be applied at every acceleration time 
throughout the maneuver, is calculated using the equation 

A Verr magbias = (bias mag )(1 0 sec )R V gailss (11) 

where bias mag is 1.5 pg, or 1.5 * 9.81 * 10’ 6 m/s 2 , and R V gailss is a random number value 
generated from a Gaussian distribution of mean equal to 0 and standard deviation equal to 1 . 

The error due to a scale factor on the delta-V is shown in the equation 

A Verr mag _ ScaleFactor (0 = ScaleFactor Var0n ||d(/)||(l 0 sec)i? V gauss (12) 

where a(t ) is the true average acceleration corresponding to time t in the acceleration file and the 
scale factor value to be used at each acceleration time is varied in a sinusoidal fashion: 

ScaleFactor Var0n = ScaleFactor + (3 xlO' 6 )(o.5°C)sin(^ L At) (13) 

where ScaleFactor is 0.01625%, T is the period of the temperature sinusoidal variation which is 
60 minutes for this analysis, and A t is the time that has elapsed since the maneuver start time. 

The magnitude error due to temperature variations assuming a 3 pg/°C bias sensitivity is 

A V6fy mag T emp y ar (0 = (3 X 10“ 6 X 9 - 8 1 mis 2 )(0.5 o C)sin(^ A t) (14) 
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The drift error due to random walk is computed at every acceleration time within the 
maneuver: 

A Verr dnft (0 = A Verr drift (tprev ) + { dri ftmag V gauss (0 (15) 

where t is the previous acceleration time, A Verr dHfi (t prev ) is the drift error from the previous 
acceleration time within the maneuver, and RV gauss (t) is a random number value generated at 
every acceleration time in the maneuver. For this analysis, t - t is 10 seconds. 

The delta-v direction bias that is applied at every acceleration time throughout the maneuver is 
determined using the cone-and-clock method. The angle of the cone, i.e. the angle between the 
acceleration vector and the error vector as shown in Figure 15 , is 

0 = cone angle = \ [dir mag )r V gauss | ( 16 ) 

The clock angle, i.e. the point on the circle formed by the cone that will be selected as the 
direction of the error vector as shown in Figure 15 , is given by 

(/) = clock angle = (360 °)RV uniform ( 17 ) 

where RV uniform is a random number value generated from a uniform distribution on the open 
interval 0 to 1 . 


= clock angle 
a 
a 


Figure 15. Direction Error Cone. 

The reference frame in which the error is computed must be defined in inertial space at every 
acceleration time within the maneuver. Three unit vectors are defined: 1) u x , along the direction 

of the acceleration in the inertial frame, 2) z7 2 , along the cross product between the first unit 
vector and the orbital angular momentum vector unit vector, and 3) u 3 , along the cross product 

between the first two unit vectors. Thus the unit vector of the acceleration vector with direction 
error at a particular time is expressed in the inertial frame as 

«„_^«.,=\a,u 1 u,]R l (4,)R 1 (0)[ 1 0 0] T (18) 
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where R l ((/>) is a rotation matrix about the first axis, R 2 {0 ) is a rotation matrix about the second 
axis, and \u x u 2 ii 3 ] is the rotation matrix that transforms the unit vector from the acceleration 
frame to the inertial frame. 

Variations Performed 

The parameters of the acceleration knowledge error model that were varied included Direction 
Bias for Formation Maintenance (FM), Direction Bias for Apogee Raising (AR), Magnitude Bias, 
Magnitude Drift due to random walk, Magnitude Scale Factor Error, and Magnitude Error due to 
Temperature Variations. Starting with the same Phase 2b 30km scenario that was used for the 
USO off nominal analysis, the thrust acceleration knowledge error model values used for the 
Baseline were 10% of the l-o values listed in Table 6. The 10% percentage was chosen because 
it produces results similar to the Phase 2b 30km baseline scenario that was used for the USO off 
nominal analysis. In each variation performed, the individual parameter value was increased by 
an order of magnitude to the l-o value, as described in Table 7. These are the actual values used 
for each parameter, i.e. they are not sampled randomly. 1 -g values were chosen because the 
difference in results between the cases were more discemable than when smaller values were 
chosen for the variations. 


Table 7. Baseline and Successive Acceleration Knowledge Error Model Variations. 


Case 

Direction 
Bias (FM) 

Direction 
Bias (AR) 

Magnitude 

Bias 

Magnitude 

Drift 

Scale 

Factor 

Temperature 

Baseline 

0.0681° 

0.0447° 

0.15 Mg 

2.25 gm/s/s 1 ^ 

0.001625% 

0.3 Mg/°C 

Direction Bias 
(FM) increased 

0.681° 

0.0447° 

0.15Mg 

2.25 gm/s/s 1/2 

0.001625% 

0.3 Mg/°C 

Direction Bias 
(AR) increased 

0.0681° 

0.447° 

0.15 ms 

2.25 gm/s/s 1 ' 2 

0.001625% 

0.3 Mg/°C 

Magnitude Bias 
increased 

0.0681° 

0.0447° 

1.5 Mg 

2.25 Mm/s/s 1/2 

0.001625% 

0.3 Mg/°C 

Magnitude Drift 
increased 

0.0681° 

0.0447° 

0.15 ms 

22.5 nm/s/s 1/2 

0.001625% 

0.3 Mg/°C 

Scale Factor 
increased 

0.0681° 

0.0447° 

0.15 Mg 

2.25 gm/s/s 1/2 

0.01625% 

0.3 Mg/°C 

Temperature 

increased 

0.0681° 

0.0447° 

0.15 Mg 

2.25 gm/s/s 1/2 

0.001625% 

3 Mg/°C 

All values 
increased 

0.681° 

0.447° 

1.5 Mg 

22.5 pm/s/s 12 

0.01625% 

3 Mg/°C 


Phase 2b 30km Formation Maintenance Maneuver Results 

All results below are for spacecraft MMS 1, which is representative of the sensitivity trends 
seen for the other spacecraft. The clock error, RSS position error, and RSS velocity error for the 
baseline case are shown in Figures 16 through 18. 
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Figure 16. Baseline Clock Bias, Drift, Acceleration Error and Covariance. 
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Figure 18. Baseline RSS Velocity Error and Covariance. 

Figure 19 shows the maximum RSS position error for all Phase 2b 30km cases described in 
Table 7, with the same trends seen in maximum RSS velocity error for each case. These values 
reveal that the position and velocity error are most sensitive to the magnitude drift. 
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Figure 19. Maximum RSS Position Error for Phase 2b 30km Cases. 

The increase of the magnitude drift parameter results in very slightly modified clock error, 
indicating a minimal sensitivity of the clock error to acceleration knowledge errors. Figures 20 
and 21 show the RSS position and velocity errors with the direction bias parameter increase as 
defined in Table 7. The increase of the magnitude drift parameter had a significant impact on the 
position and velocity errors over the propagation period immediately following the second 
maneuver. The position and velocity error was significantly different only in the time between 
the second maneuver and the following perigee passage, and to a much smaller degree the orbit 
following this perigee passage. 
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Figure 20. Position Error and Covariance with Magnitude Drift Increased to l-o. 
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Figure 21. Velocity Error and Covariance with Magnitude Drift Increased to l-o. 


Phase 2b 400km Formation Maintenance Maneuver Results 

Phase 2b 400km was considered in addition to Phase 2b 30km because the maneuvers were 
significantly larger than those in Phase 2b 30km, which allowed for an investigation of how the 
size of the maneuver affects the sensitivities. The size of the maneuvers for the Phase 2b 400km 
scenario were on the order of 6 to 11 meters per second, versus those for the Phase 2b 30km 
scenario, which were on the order of 1 to 3 meters per second. The same sensitivity patterns were 
seen for Phase 2b 400km as were seen in Phase 2b 30km, except that the direction bias increase 
did not produce a significant change in the maximum RSS position or velocity error values. 

Phase 2a Apogee Raising Maneuver Results 

The thrust acceleration knowledge error for the longer Apogee Raising maneuvers in Phase 2a 
were also varied using the same method as for Phase 2b 30km described previously, however it 
was determined that due to the fact that the Apogee Raising maneuvers occur when the spacecraft 
is receiving GPS measurements, the impact of thrust acceleration knowledge errors is minimal. 
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CONCLUSION 

This study indicated that increases in the clock drift random walk parameter h_ 2 produces the 
largest increases in the navigation error among the USO clock error model parameters, with a 
secondary sensitivity to the amplitude of the umbra temperature events b 3 . The maximum clock 
bias error observed in each primary clock error model variation was well below the maximum 
clock bias error requirement of 325 microseconds, indicating significant margin. For the thrust 
acceleration knowledge error model parameters, it was determined that for both shorter and 
longer maneuvers, the total navigation error is primarily sensitive to knowledge magnitude drift 
due to velocity random walk during the maneuver. However, this sensitivity is only significant 
when there are long periods of propagation following the maneuver. 
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